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Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann 
transport terms in the distribution function equations for electrons and holes, we derived a closed 
set of diffusion equations for carrier densities and temperatures with self- consistent coupling to 
Maxwell’s equation and to an effective optical polarization equation. The coherent many-body effects 
are included within the screened Hartree-Fock approximation, while scatterings are treated within 
the second Born approximation including both the in- and out-scatterings. Microscopic expressions 
for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the 
momentum and energy relaxation rates. These rates expressed as functions of temperatures and 
densities lead to microscopic expressions for self- and mutual- diffusion coefficients in the coupled 
density-temperature diffusion equations. Approximations for reducing the general two-component 
description of the electron-hole plasma (EHP ) to a single-component one are discussed. In particular, 
we show that a special single-component reduction is possible when e-h scattering dominates over 
c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show 
that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion 
coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions 
lead to new perspectives into the roles played in the single-component reduction by the electron- 
hole correlation in momentum space induced by scatterings and the electron-hole correlation in real 
space via internal static electrical field. Finally, the theory is completed by coupling the diffusion 
equations to the lattice temperature equation and to the effective optical polarization which in turn 
couples to the laser field. 
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I. INTRODUCTION 


A typical semiconductor-based optoelectronic device, 
such as a diode laser, can be considered as consisting 
of three subsystems: an optical field, an electron-hole 
plasma (EHP), and a host crystal lattice. Light gen- 
eration, propagation, amplification, and diffraction de- 
termine the behavior of the optical field, while electri- 
cal conduction, plasma diffusion, and carrier generation 
and recombination determine that of the plasma. The 
host crystal is very often represented by various phonon 
modes, with the longitudinal optical (LO) phonon mode 
being the most important one for optics and carrier trans- 
port of III-V semiconductors at room temperature. Obvi- 
ously the whole system involves the interplay of optical, 
electrical and thermal processes. From the perspective 
of conversion and conservation, energy is stored in the 
form of photon energy, kinetic and thermal energy of the 
EHP, and thermal energy of the host crystal represented 
by the phonons excited and annihilated. Though various 
processes occur in different space and time scales, their 
couplings should be treated in a self-consistent fashion 
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to arrive at correct coupled equations, or to appropri- 
ately decouple the equations in various limiting cases. 
While modeling and simulation of semiconductor lasers 
and other optoelectronic devices typically use the rate 
equations [1] or the semiconductor Bloch equations [2] 
that treat only the coupled electrical-optical subsystems, 
other approaches have focused on the thermal aspect [3]. 
In Refs. [4, 5], a self-consistent approach was attempted 
to combine all three subsystems and the relevant pro- 
cesses based on earlier work on plasma heating [6-9] in 
semiconductors. The approach is, however, valid only for 
a single-mode laser, or a spatially homogeneous laser. 

There is plenty of manifestation of spatial inhomo- 
geneities in a semiconductor device. Laser beam fila- 
mentation, dynamic beam steering, and multiple trans- 
verse mode formation and competition are some of the 
important examples where inclusion of spatial, or spatio- 
temporal variation is necessary. Additionally, with the 
inclusion of temperature variables, more spatial phe- 
nomena can be described, such as thermal lensing and 
formation of hot spots, and catastrophic optical dam- 
age (COD) in diode lasers. For advanced semiconduc- 
tor lasers such as vertical-cavity surface-emitting lasers 
(VCSELs), or Master-Oscillator coupled with Power- 
Amplifier (MOPA), spatial inhomogeneity is a more 
prominent issue when efficient coupling of multimode 
VCSELs to multimode fibers are desired, or when care- 
ful transverse mode engineering is required for certain 
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applications. In terms of theoretical description, spatial 
inhomogeneity of a semiconductor laser is usually dealt 
with using Maxwell’s equation with diffraction terms, 
coupled with carrier diffusion equation [10]. This set of 
equations, which is a direct generalization of the rate 
equations to the spatially inhomogeneous case, can be 
further augmented by adding the material polarization 
equation to account for gain and refractive index dis- 
persion [11]. This approach has recently been used to 
simulate MOPAs [12, 13] and VCSELs [14]. An alter- 
native approach is to start from the Boltzmann- type of 
transport equation for the carrier distribution functions 
as described in Ref. [15]. This approach is a very general 
one, but it involves quite extensive computation, since 
A>resolved interband polarization equations have to be 
solved self-consistently for all space and time grid points. 
For nearly all practical purposes, lasers and other opto- 
electronic devices operating on time scales longer than 
a picosecond can be more efficiently dealt with by the 
moment-equation type of approach. To account for both 
carrier density and carrier energy (thus carrier temper- 
ature) inhomogeneities within the framework of the mo- 
ment equations, one needs to derive a set of partial differ- 
ential equations for up to the second order moments from 
the corresponding Boltzmann transport equations. Such 
an approach, which is alternatively called hydrodynamic 
approach, will lead to a generalization of the single-mode 
laser model in Ref. [4] to spatially inhomogeneous cases 
including thermal and heating effects. In the past, car- 
rier diffusion equation has been derived from the Boltz- 
mann .transport equation in combination with Maxwell’s 
equation and the optical polarization equations by sev- 
eral authors [15, 16]. Since only the zeroth-order mo- 
ment equation is derived [15], plasma heating cannot be 
described. This lack of systematic derivation of the tem- 
perature equation creates problem in correctly choosing 
the energy or temperature diffusion coefficient [17]. Fur- 
thermore, whenever two or more types of diffusion pro- 
cesses exist, mutual diffusions, or cross diffusions, occur. 
The quantitative significance of these processes needs to 
be examined, which is only possible with a systematic 
derivation of the coupled diffusion equations and all the 
corresponding diffusion coefficients (DCs). 

This work sets out to derive such a set of coupled 
macroscopic equations for carrier densities and carrier en- 
ergies from the coupled Boltzmann-Bloch transport equa- 
tions using moment-equation approach. We pay special 
attention to the treatment of carrier-carrier (c-c) scat- 
tering and carrier-LO phonon (c-LO) scattering. Dif- 
ferent from previous approaches in this regard, micro- 
scopic electron-hole (e-h) and c-LO scattering terms in 
the Boltzmann transport equations are used directly to 
obtain the corresponding energy and momentum relax- 
ation rates. Various DCs are then expressed in terms of 
these momentum relaxation rates. The energy relaxation 
rates are used to describe energy exchange between dif- 
ferent subsystems, which translate to temperature vari- 
ations consequently. Such an approach allows a detailed 


study of the DCs for a given quantum well structure. De- 
tailed numerical results and analysis of these coefficients 
will be presented in a subsequent article. 

The article is organized as follows. In Section II, we 
introduce the basic physical considerations and starting 
equations. This is followed by a general derivation of the 
moment equations and the treatment of the cutoff issue, 
which is central to any moment-equation based approach. 
The general form of the hydrodynamic equations for the 
general two-component situation is derived in the fourth 
section. The fifth section deals with the specialization to 
the cases of single-component- approximation where we 
will discuss the consequence of the ultrafast e-h scatter- 
ing for the drift momenta of electrons and holes. We will 
also discuss the well-known ambipolar diffusion approxi- 
mation (ADA). In the sixth section, the very important 
issue of optical polarization is treated. We use two dif- 
ferent approaches to close the hierarchy related to the k- 
resolved polarization equations. In Section VII the corre- 
sponding Maxwell’s equation and the lattice temperature 
equation are introduced to form the complete set of equa- 
tions for a more complete description of semiconductor 
lasers. Though we use laser as an exemplary device in 
this work, the treatment and the resulting equations can 
be used for other optoelectronic devices, such as photo- 
conductors [18] and photodetectors [19], with little or no 
modification. In Section VIII we discuss and comment on 
some general aspects of our theory and sum up the major 
results of our article before we present detailed consid- 
erations in the appendices for scattering rates and DCs. 
Appendix A deals with momentum and energy relaxation 
rates due to e-h scattering, while the corresponding rates 
due to c-LO phonon scattering are treated in Appendix 
B. In Appendix C, we list DCs for the two-component 
case. Finally, expressions for all the DCs under the one- 
component approximation are given in the last Appendix. 


II. MODEL AND BASIC EQUATIONS 

We begin this section with some general remarks about 
the spirit of our approach. As is often assumed, laser- 
matter interaction in a semiconductor laser is described 
with an EHP model. The standard argument to sup- 
port such a plasma model is that the typical room- 
temperature lasing density is around 10 12 cnU 2 , well 
above the Mott density. W T hile we use this model in this 
article, its adequacy in a spatially inhomogeneous sys- 
tem needs special scrutiny. We comment on this in the 
final section of this article. Within the plasma model, 
Coulomb interaction, is usually characterized by an ex- 
citonic enhancement of the optical transition and car- 
rier density-dependent bandgap renormalization [2, 16] in 
the coherent part of the semiconductor Bloch equations 
(SBEs). In addition, this relatively high density leads to 
ultrafast carrier-carrier scattering within 100 femtosec- 
onds, which dominates carrier dynamics on the short 
time scale and affect the interaction of the EHP with 
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a laser field in several important ways: First, these ul- 
trafast collisions thermalize carriers in properly popu- 
lated subbands [20], which justifies the standard assump- 
tion of quasiequilibrium for carrier distributions when 
longer time scale dynamics is of interest. Second, as will 
also be shown later in this article, ultrafast e-h scatter- 
ing correlates populated subbands such that the whole 
EHP can be characterized with a single temperature and 
drift velocity, but different chemical potentials for dif- 
ferent subbands under normal conditions. This under- 
standing will be elaborated in the subsequent sections 
and in appendix A by examining in detail the effects 
of the scattering on momentum and energy relaxations. 
Due to these ultrafast collision events, local quasiequi- 
librium is established in less than a picosecond. This 
allows carrier distribution functions to be described lo- 
cally by Fermi-Dirac distributions with finite drift mo- 
menta, while spatial inhomogeneity is taken into account 
by the space-dependent macroscopic “parameters” (den- 
sities and temperatures) of such distribution functions. 
Such a treatment further allows other slower processes 
be incorporated in a hydrodynamic fashion. Reasoning 
and understanding of this type underpin the moment- 
equation approach adopted in this work. 

Specifically, we consider an intrinsic semiconductor 
quantum well of width w in the z direction and of area 


S in the x-y plane. The inhomogeneity occurs in the 
plane of the quantum well layer, while fixed profiles for 
carrier distribution and optical modes are assumed in 
the vertical direction to the quantum well due to tight 
quantum confinement and optical wave-guiding. The ex- 
tension to include the vertical inhomogeneity can be nat- 
urally made when we include a vertical transport model, 
such as the one used in Ref. [21]. This is beyond the 
scope of the present article. In a typical edge emitting 
laser (EEL) , the inhomogeneity appears in the plane of 
the light propagation. For VCSELs, this plane is trans- 
verse to the light propagation. Our starting point for 
the semiconductor medium is the semiconductor Bloch 
equations [16, 22], generalized for the spatially inhomo- 
geneous case in Ref. [15, 23], which may now be called 
Boltzmann-Bloch equations. The possible local charge 
unbalance requires that the Poisson equation be included. 
The complete set of equations required for such a system 
is therefore the Maxwell-Boltzmann-Bloch-Poisson equa- 
tions (MBBP) [2, 15, 16, 22-24] for the non-equilibrium 
distribution functions n a (k , r) (a = e, h for electrons and 
holes, respectively), interband polarization p(fc,r), elec- 
tric potential 4>(r), and the laser field £(r,£), with k and 
r being the 2D vectors in k space and real space, respec- 
tively. The MBBP equations are collected as follows: 
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$dl£{r,t)-V£{r,t) = -^dnv a + V b ) , ( 1 ) 

d t n a (k, r ) + \d k e°{k,r) - d r n a {k,r) - \d r r) + «°#(r)] • d k n°{k, r) 

= R a (k,r) + d t n a (k,r)\ co t , . ( 2 ) 

d t p{k,r ) = -{ [£ e (k,r) + e h (-k,r)]p(k,r)-iQ(k,r) [n e {k,r) + n h (-k,r) - l] + d t p(k, r)| coi , (3) 

d 2 Mr) = -^E a , k 9 a ^,r), ( 4 ) 


with various terms given below: 
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R a (k,r) = A a (fc,r) -7„ r n Q (fc,r) - B sp n e {k, r)n h (k,r) + i [fi(fc,r)p*(fc, r) - n*(fc,r)p(fc,r)] , 

e a (k,r ) = e a (k) + q°$(r) + Se a (k,r ) , 

6e a {k,r) = - r)V Sik - k ' +<5 a ,v “ v k-k’) , 

hQ(k,r) = n(k)£(r,t) + J2 k ' P{^' > r Ws,k-k' , 

B a — y Sfc [a (k)p(k, r) + c.c .] , 
d t n Q {k,r)\ co i = d t n a (k,r)\ aa + d t n a (k,r)\ eh + d t n°(k,r) \ LO , 
d t p{k,r)\ co i = — r o (fc)p(fc, r ) + £ fc , r*(fc, k')p(k',r) , 


(5) 

( 6 ) 

(7) 

(8) 
(9) 

( 10 ) 

( 11 ) 


where d t = d/dt, d T = d/dr , d k = d/dk, d\ = d 2 /dt 2 , is the renormalized carrier energy, 5 e a (k,r ) is the cor- 

y2 = Q2 + Q2 j dz 2 ) Q2 = q 2 j dx 2 + q 2 j Q y 2 A i S0) £ a r ) rec tion to the single-particle carrier energy e°(k) due to 
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exchange interaction and due to the Coulomb-hole self- 
energy. £l(k,r) is the renormalized Rabi frequency. V $i k 
is the screened Coulomb potential, for which we use the 
single plasmon pole model in this work, and Vk is the bare 
one. In addition, fi(k) is the interband optical dipole ma- 
trix element between electron state \ch) and hole state 
|x? —k). Vb is the optical polarization of the unexcited 
semiconductor, while V a accounts for the electronic con- 
tribution from photoexcitation, so the total material po- 
larization V = V a +V b . Their treatment together with 
Maxwell’s equation is detailed in Section VI and Sec- 
tion VII. Furthermore, h is the Planck constant, q a = qpe 
is the carrier charge for electrons and holes, £q is the elec- 
tric constant, is the relative permittivity of the unex- 
cited semiconductor, V = wS is the volume of the ac- 
tive region, S QiV is the Kronecker delta, and 2 J2 k means 
summation over all allowable momentum (fc) states, in- 
cluding spin (account for the 2), for a single subband. 
Additionally, the lumped generation-recombination (g-r) 
contributions R a (k,r) consists of current injection term 
A Q (fc,r), non-radiative recombination term with coeffi- 
cient 7 nr , spontaneous emissions term with coefficient 
B sp , and stimulated interaction term involving the renor- 
malized Rabi frequency f 2(fc, r). Detailed treatment of g- 
r term is available in Section VI. Finally, collisions! con- 
tributions denoted by subscript col lead to decay in the 
interband polarization dtp{k,r)\ co i and relaxation in the 
carrier distributions dtn® (fc, r)| co j. We include explic- 
itly both the so-called out-scattering term r o (fc) and the 
(non-diagonal) in-scattering term F, (/c, /c'). The domi- 
nant scatterings considered in this work are c-LO, a-a, 
and e-h scatterings, as separately denoted in Eq. (10). 

Before w r e start the formal derivation of the moment 
equations, a few r remarks are in order: First, parabolic 
band is assumed for electrons and holes with effective 
mass m a , respectively. This approximation is valid wffien 
the well width (u>) is small enough and plasma density 
is not too high, such that higher subbands are not pop- 
ulated and thus ignored. This means that carrier ki- 
netic energy, given as e a (k) = h 2 k 2 /2m ai is proportional 
to the second order moment. As w^e will see later, this 
simplification allows direct association of the second or- 
der moment with the total carrier kinetic energy and 
the first order current with this energy partially. As 
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such, the extension of the following derivations to the 
case of non-parabolic bands is not straightforward. Sec- 
ond, all the scattering terms will be treated explicitly in 
the moment equations, without resorting to relaxation 
rate approximations or leaving them at the formal level. 
This is where our approach differs from those earlier ap- 
proaches [15, 16] in deriving the moment equations; this 
will become more evident as -we proceed with the deriva- 
tion. Finally, short-hand notation will be adopted for 
convenience and brevity, such that n k stands for n a (k , r), 
Se Q for Se a (k,r ), $ for $(r), and so on, in the remaining 
part of this article unless indicated otherwise. 

HI. MOMENT EQUATIONS AND CUTOFF 

We begin with introduction of the moment and current 
of the n-th order associated with the non-equilibrium dis- 
tribution function n k as follows: 


C( r ) = 

> 

* fc 

(12a) 

Jn(r) = 

* fc 

(12b) 


wffiere F®, denoting the n-th order weight function, are 
1, hk, and h 2 k 2 /2m a for n - 0,1,2, respectively, and 
v k = 9k€ k /h — hk/m a for parabolic bands. We con- 
sider up to the second order for moments. Note that 
only the trace of the second order moment tensor needs 
consideration here, and it relates to energy. For clar- 
ity, we represent the first three considered moments and 
currents with conventional symbols: density N a = 
momentum P a = energy E a = t/£; density cur- 
rent = J 0 a , momentum current = J®, energy 
current J% = . Jp is a tensor despite its mislead- 

ing vectorial notation. As is customary (see, e.g. Ref. 
[25]), we derive the moment equations by summing over 
all degrees of freedom, t.e., applying (2 /S)Yj k on fhe 
Boltzmann transport equation [Eq. (2)] with the corre- 
sponding weight function F®. It is straightforward to 
show that the first three moment equations for electrons 
and holes can be w r ritten as belowx 


d t N a + d r ■ = R% , (13) 

d t P a + d r Jp+N a dr(5t a + q a $) = Rp + d t P a \ eh + d t P a \ LO , (14) 

d t E a + dr ■ J% + d r (<5e“ + </"$) ■ J%- = R% + d t E a \ eh + d t E a \ LO , (15) 


where we have neglected the weak Jfc-dispersion in the the corresponding weight functions. For example, 
many-body correction 5e° . Terms on the right-hand side 
above are the result of summing the corresponding terms 
in Eq. (2) over all degrees of freedom after multiplied by 
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Rn,P,E — C ^ / FN, p ,ER°‘(k’ r ) ■ (16) 

k 

Intuitively, a-a scattering does not change the total car- 
rier number, momentum, and energy for each band, so 
they vanish from the above equations. Furthermore, e-h 
and c-LO scatterings do not alter the total carrier num- 
ber within each band, but they do survive the summation 
and remain in Eqs. (14-15). 

While the formal derivation of the above moment equa- 
tions [Eqs. (13-15)] is exact and straightforward, these 
equations are not in a closed form yet. There are sev- 
eral 'reasons for this: (1) the first order current J% is not 
completely given by the first and second order moments 
(pa and E a ) in the presence of anisotropy, where equa- 
tions for the other (trace-less diagonal and off-diagonal) 
elements of the second order moment are required; (2) 
the second order current is connected to higher order 
moments; (3) the terms on the right-hand side depend on 
summation over the carrier distribution functions, thus 
are not fully known in terms of the first three moments 
and currents. These are well known reasons that exist 
in general for the moment-equation approach, no matter 
whether it is applied for transport problems in microelec- 
tronics or in fluid dynamics, which lead to the so-called 
hierarchy problem. In addition, there is an extra com- 
plication in optical problems as we are considering here: 
The carrier distribution functions are coupled to the in- 
terband polarization p(fc,r), another space and carrier 
momentum-dependent distribution function. As a result, 
we are essentially dealing with coupled Boltzmann trans- 
port equations for three distribution functions: n e (fc,r), 
n h (k,r), and p(fc,r), though all transport terms involv- 
ing explicit'spatial variation of p(fc, r) are ignored [23]. In 
general, moment equations for all these distribution func- 
tions should be sought. To obtain a closed form of equa- 
tions for these macroscopic quantities, two approaches 
are typically used: One is to derive the lower order of the 
moment equations and cut off the hierarchy by setting the 
higher order moments to zero. As it is evident, the cou- 
pling to the p(k, r ) 1 s makes such moment cutoff impracti- 
cal. Another approach is to assume that the distribution 
functions are well approximated by known distributions 
characterized by some macroscopic parameters. These 
macroscopic parameters can be linked to those moment 
variables. In fluid dynamics or in microelectronics, for ex- 
ample, one assumes that the system is locally described 
by the drifted Maxwell distribution [26]. Similarly, in 
semiconductor laser theory, the quasiequilibrium condi- 
tion is well established [22, 24, 27]. Quasi-equilibrium 
here means that electrons and holes, driven out of mu- 
tual equilibrium by laser field and external pumping, are 
separately characterized by the equilibrium distributions 
of each subsystem in their inertia frame of reference. The 
physical mechanism responsible for the establishment of 
this quasiequilibrium is the ultrafast a-a scattering on 
the femtosecond time scale [28]. Using this quasiequi- 


librium assumption and neglecting the other elements in 
the second order moment tensor except its trace (valid for 
isotropic physical systems), we can readily truncate the 
hierarchy associated with n Q (fc,r). We will come back 
to the hierarchy problem associated with p(fc,r) later. 

First of all, we assume that the quasiequilibrium distri- 
butions of the EHP are given by the drifted Fermi-Dirac 
(DFD) distribution functions, 

K = fk-k% = I 1 + ex P [&>( e fc-fcS - Vf)] } , (17) 

where fcg is the drift wavevector and pf is the chem- 
ical potential. Moreover, = l//c#T a , where ks is 
the Boltzmann constant and T a is the temperature for 
electrons or holes. The drift wavevector is related to 
the first order moment and pf is given by (3 a pp = 
In [exp (7r/3 Q ^ 2 A TQc /m Q ) - l] in 2D case. We note that 
three parameters are needed to characterize a DFD func- 
tion given by Eq. (17). A total of six parameters for the 
electron and hole distribution functions can be uniquely 
associated with the six moment variables. With the aid of 
the known functional form of the DFDs, the right-hand- 
side terms in Eqs. (13-15) are calculated as functions 
of carrier densities, drift wave vectors, and temperatures. 
Finally, using the definition of moments and currents, 
Eqs. (12), with n£ replaced by the DFD function, we can 
show that currents of the first three orders depend on the 
moments as follows: 

= N«u a , (18) 

J«=u a P a +W a I > (19) 

J% = 2u a W Q + ±u a u a • P Q , (20) 

where I is the unit tensor, 

< 2I > 

k 

is the thermal part of the carrier energy E a , and 

E a = W a + • F° ■ (22) 

Above, P a = N a hkp and a drift velocity is introduced: 


u° = P a /m a N a = hkp/m a . 

Obviously the carrier energy E a is related to temperature 
via its thermal energy part given by Eq. (21). Therefore, 
these relations between currents and moments connect all 
dynamical variables in a closed form, leading to a closed 
set of equations for {N a , P a , W°}. 
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IV. GENERAL TWO-COMPONENT MOMENT between the moments and currents given in Eqs. (18-22). 

EQUATIONS The resultant equations are written more specifically as 

follows: 

The moment equations given by Eqs. (13-15) are now 
closed for {N a , P a , W a } after applying the relations 

I 


d t N a + d r • ( u a N a ) = R% , ( 23 ) 

d t P a + d r ■ ( u a P a ) + d r W a + N a dr {6e a + 9 “$) = R a p + d t P a \ eh + d t P a \ LO , (24) 

d t W Q + d r ■ (2 u a W a ) - u a ■ d r W a = + d t W a \ eh + d t W a \ LO , (25) 


where 2 u a W a = Jyy is the thermal energy current den- 
sity. In addition, 

Rfr =R%+ \m a u a • u a R% (26) 

dtW*\ eh = dtE* | efc - ll" • d t P a \eh , (27) 

- 3t£ Q |io - u* • • (28) 

Equations (23-25) are the general form of the moment 
equations which describe the lateral spatial-temporal car- 
rier dynamics in a semiconductor laser. These equations 
should be solved together with the Poisson equation, 
Eq. (4), the still fc-resolved polarization equations, and 
Maxwell’s equation. The scattering terms in the above 
equations are specified in their general form or in the lin- 
earized form in Appendices A and B. From Eqs. (26-28), 
we see that terms {Rfo, d t W a \ e h, dtW Q \z,o} differ from 
{R%, d t E a \ ehy d t E a \io} correspondingly by nonlinear 
terms in P a . As it turns out (see Appendices A and B), 
the nonlinear terms drop out in the thermal energy re- 
laxation expressions. Finally, it can be shown that these 
terms vanish for Rf v as well if translational invariance in 
k space is assumed for the lumped g-r term R a (k,r ) in 
Eq. (2), which intuitively makes sense as internal energy 
generation and recombination should be independent of 
translational momenta. 

To simplify the above set of equation further, we in- 
tend to eliminate the momentum equation (24). First, we 
point out that the g-r term Rp is negligible as compared 
to the two scattering terms since momentum relaxation 
is dominated by ultrafast scattering events, given that 
no appreciable momentum transfer occurs accompany- 
ing the generation and recombination processes. Fur- 
thermore, as shown in Appendices A and B, the two 
scattering terms can be linearized, as a consequence of 
the assumption of quasiequilibrium for the electron-hole 
plasma. Possible violation of this assumption comes from 
the presence of a strong electric field in the quantum well 
plane that tends to drive the system out of quasiequi- 
librium and into a nonlinear regime. Substituting the 
scattering terms by their linearized forms as given in 
Eqs. (A8) and (B6), the momentum equation (24) is now 
rewritten as follows: 


d t P a + d r (u a P«) + d r W a +N a d r {6e a + q Q $) 

= -m r j eh (tt“ - u°) - 7 l 0 P a , (29) 

where a ^ ft, rn r y e h/Tn Q N Q corresponds to the mo- 
mentum relaxation rate due to e-h scattering, m r = 
77t e 77ih / (7n e + rrih) is the reduced mass, and Jiq is the 
momentum relaxation rate due to c-LO scattering. To 
obtain an explicit expression for the density current for 
each carrier type, we adopt the adiabatic elimination ap- 
proximation for the above equation [16]. For a weakly 
inhomogeneous system, we can further ignore the nonlin- 
ear (second on the left side) term in P a . The resulting 
solution for the momentum is given as 


P a = —m a n a Y, [8rW v + ^ 3 , {Se p + g !/ 4>)] 

i/=e,/i 

—m a fi a T] a [d r w a + N a dr (Se a + q a *)] , (30) 
where the two factors are defined by 


Hot 


Vc* 


'Yeh 


7 LO^Lo( m e + m h ) + % h Ep=e t h m ^LO 


* 31 ) 


7 LO m e ± m h 

7 


(32) 


where a,/? € {e, /i) | a ^ ■/? and y^ h = 7 e /i/A r ^. The 
two terms in the above equation have distinct physical 
meanings: The first term relates to the e-h scattering 
that tends to equilibrate the two carrier types, while the 
second term describes the equilibration process between 
each carrier type and the LO phonon subsystem. We 
mention that similar formulation and results for the e-h 
scattering have been derived previously in the study of 
the negative mobility of minority carriers in semiconduc- 
tor quantum wells [29-34]. 

Finally, equations for temperatures are very often pre- 
ferred over those for thermal energies. Fortunately, a 
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unique transformation exists under the DFD approxima- 
tion, since W a = W a [fi a F (N a ,T a ),T a ] = W a (N a ,T“). 

It is given below: 

d t T a = fadtW* - j%d t N a , (33) 

and the resultant temperature equation goes as: 

d t T a + dr • J? - ■ drW a + drj% • ~ d r J?y ■ 

= 3w (Rw + 9tW a \ LO + d t W a \ eh ) - j%R% , (34) 
where the temperature current Jf = jw^w ~3 n^n an< 3 

jfr = (^T-^lAT-r 1 , (35a) 

j% = d N °W a \ T ° (3r-W’°Vr 1 • (35b) 

To summarize this section, Eqs. (23,34) form the closed 
set of diffusion equations in terms of carrier densities 
and temperatures for each component after we replace 
momentum P a or u a in Eqs. (23,34) with its adiabatic 
solution of Eq. (30). The corresponding density and tem- 
perature currents, and are given by 

J%r = P a /m a , (36) 

j'i- = [2 jfr ( W/N * ) - m J a N • (37) 

The above currents contain the gradients of four macro- 
scopic variables: __N a and T a for a = e,/i. In general, 
we can introduce a 4 x 4 diffusion matrix, Dxy ? with 
X, Y e {N e , N h ,T e >T h }. While the diagonal elements 
of this matrix represent the self-diffusion coefficients, the 
off-diagonal elements denote various mutual- or cross- 
diffusion coefficients. A complete list of all the coeffi- 
cients is given in Appendix C. 

V. SINGLE-COMPONENT APPROXIMATION 

In this section, we consider two approximations that al- 
low reduction of the general two-component description 
of the electron-hole plasma to a single-component one: in 
the limit of strong electron-hole scattering and the often 
used ambipolar diffusion approximation. Detailed analy- 
sis is provided for the comparison of the two approaches 
and they are found to produce consistent results in the 
linear regime of description. 

A. The Limiting Case of Strong Electron-Hole 
Scattering 

While we used the quasiequilibrium assumption in the 
above derivation, which is due to ultrafast carrier-carrier 
(a- a) scattering within each band, e-h scattering, which 
is on the same time scale [29, 35], has been retained 


in Eqs. (13-15), together with c-LO scattering (in sub- 
picosecond range) and carrier diffusion and energy trans- 
port (in nanosecond range). Thus self-consistency de- 
mands that we further consider the dynamical correlation 
between electrons and holes imposed by e-h scattering. 
In this subsection, we take on the issue of how e-h scat- 
tering reduces the general two-component description to 
a single-component one for the EHP near quasiequilib- 
rium. As shown in Appendix A, detailed balance (DB) 
requirement for quasiequilibrium in the sole presence of 
e-h scattering leads to these conditions: 

T e - T h , (38) 

= ti* , (39) 

which is intuitively apparent as electron-hole scattering 

is meant for equilibration between the two carrier types. 
The above conditions are the same as the a posteriori 
requirements needed for the ambipolar diffusion approx- 
imation [36], and they now permit us to settle the is- 
sue of reducing the original two-component problem to 
a single-component one if the EHP is initially neutral in 
real space. Further discussions along the line of standard 
ADA reduction will be presented in the next subsection. 

To examine this issue in greater detail, we consider the 
dynamics around the DB state by looking at the equa- 
tions for momenta and energies with scattering terms 
linearized around the DB state. The linearized scat- 
terings terms are derived in Appendix A and the cor- 
responding momentum equations are given by Eq. (29). 
As we see from the adiabatic solution to these equations 
[see Eq. (30)], in general, Eq. (39) is not valid. This 
means that the DB is not sustainable and the correspond- 
ing single-component reduction does not hold. However, 
when t £ 0 7 ^, the second term in Eq. (30) can be 

neglected. The omission of the second term conveniently 
leads to the conclusion of u e = u h = u y if N e = N h = N. 
Thus 

u = -■£ [W + Nd r {Se v + <?"$)] , (40) 

i/=e,h 

where p — p e = Vh as defined in Eq. (31). Take a note 
of the difference between this line of single-component 
reduction and the standard ADA line, as presented in 
next subsection. Here the drift velocities are equal, ir- 
respective of the internal electric field (— <9 r $). In other 
words, internal field is not required to maintain equal ve- 
locities. The ultrafast e-h scattering alone maintains the 
charge neutrality if the system is neutral initially. It is 
clear from the above discussion that the validity of equal 
drift velocities for the two components of unequal masses 
requires that e-h scattering dominate carrier-LO phonon 
scattering. This is intuitively easy to understand from 
the physical point of view. The role of e-h scattering is 
to correlate electrons and holes dynamically and equili- 
brate their drift velocities. By contrast, c-LO phonon 
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scattering is to generate different individual drift mo- 
menta because of unequal masses and therefore unequal 
scattering rates with LO phonons. When the e-h scat- 
tering strength is comparable to that of c-LO scattering, 
the unequal c-LO phonon scattering rates for electrons 
and holes will be enough to counteract the homogenizing 
role played by e-h scattering, thus resulting in different 
drift velocities. At this point, one may argue if the DB is 
still a valid concept in this situation. We point out that 
the DB condition as expressed for e-h scattering alone in 
Appendix A is no longer true. Rather, the DB between 
in- and out-scatterings for a given fc- state must include c- 
LO phonon scattering a s well, as the latter becomes non- 
negligible. Since c-LO phonon scattering is sensitively 
dependent on temperature with e-h scattering being on 
density, it is clear that conditions (38-39) will no longer 
be valid for high temperature and relatively low carrier 
density. 

For now, let us continue discussion of the limiting case 
when c-LO phonon scattering is much weaker than e-h 
scattering. The hydrodynamic equality of drift velocities 
of electrons and holes means that N e = N h = N will 
be maintained if the EHP is neutral initially, according 
to Eq. (23), the continuity equation. Major g-r contribu- 
tions in R % r are the same for electrons and holes, which 
will be labeled as R jv under the single-component ap- 
proximation and examined further in the next section. 
As charge neutrality under strong e-h scattering can be 
maintained dynamically, therefore, the Poisson equation 
is automatically satisfied. The density equations for elec- 
trons and holes are reduced to a single one for the plasma 
density N: 

dtN + dr-J N = R N - (41) 

To obtain the corresponding energy equation in the 
single-component case, we notice that the temperature 
equality indicated in Eq. (38) signifies interdependence 
of the energy equations for electrons and holes. Because 
e-h scattering conserves the total energy of the EHP sub- 
system, it is natural to take the total carrier energy as 
the second dynamic variable. Applying Eq. (39) for the 
drift velocities when summing up Eq. (25) for electrons 
and holes, it is found that the total thermal energy obeys 
the following equation of dynamics: 

d t W + d r ■ J w = Rw + d t W \ LO , (42) 

where W = W e + W h , J w = 2 uW, R w = R e w + R h w , 
d t W \ LO = d t W e I LO + d t W h \ L0 . In deriving the above 
equation, we have assumed weak inhomogeneity in the 
system so that nonlinear terms have been dropped. Up to 
this point, the set of moment equations for the EHP has 
been reduced to only two: one for plasma density N and 
one for its total thermal energy W. Similar to Section IV, 
an equation for plasma temperature T — T e = T h is 
derived as follows: 


d t T + dr ■ Jt + d r ]N • Jn - drjw ' Jw 
= jw {Rw + dt\V\Lo) ~ JnRn » (43) 

where Jt — jwJw - Jn J n ■ 

Bearing great resemblance to ordinary diffusion equa- 
tions, the equations derived here, Eqs. (41) and (43), 
include many-body corrections and apply to a neutral 
EHP. To define the related diffusion coefficients, the cur- 
rents in the equations need to be expressed in terms of 
the gradients of plasma density and temperature: 

J A : = -D NN d r N-D NT d r T t (44) 

J T — —Dtn d r N — Dtt d r T * (45) 

Explicit expressions for the DCs in terms of material pa- 
rameters and thermodynamic variables, N and T, are 
given in Appendix D. 

B. The Ambipolar Diffusion Approximation 

The single-component reduction discussed in the last 
subsection is valid only when e-h scattering is much 
stronger than any other scatterings. In a typical HI-V 
semiconductor device, this is true only for the case of 
high carrier density and low temperature where e-h scat- 
tering is predominant over c-LO scattering. At room 
temperature, these two scatterings are about the same 
order of magnitude. Thus the reduction procedure above 
becomes questionable. Another single-component reduc- 
tion procedure is the so-called ambipolar diffusion ap- 
proximation. Even though it seems to us quite difficult 
to justify purely from scattering analysis the validity of 
the ADA, numerical simulation by Held et al. [37, 38] has 
indeed shown that the ADA is a quite good approxima- 
tion at nanosecond time scale and at high density. The 
difference between standard ADA and the reduction pro- 
cedure described in the last subsection can be seen in the 
following way. From Eq. (30), we notice that the total 
current consists of diffusive and conductive parts: one 
proportional to the gradient of the thermal energy and 
the other one to the field (-<9 r 4>), respectively. When 
other scatterings are negligibly small compared to e-h 
scattering, the diffusive currents for electrons and holes 
are equal when the densities are the same, so that there is 
no need for the conductive currents. In the more general 
cases, however, non-zero conductive currents are needed 
to maintain the total currents for electrons and holes to 
be the same. By requiring the electron and hole currents 
to be the same, an expression for the non-zero internal 
field is obtained. Substitution of the internal field ex- 
pression into either one of the density current expressions 
leads to the ambipolar density current: 

J a N m = ~D% m N d r N - Df$d r T , (46) 
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where the ambipolar diffusion coefficients are given be- 
low: 


s%±3. 

T^elio + m ^LO 


and 


(47) 


S« = d x W a + Nd x Se a , (X = A T , T ) . (48) 

As seen from the above expressions for the ambipolar 
DCs, they are independent of the e-h scattering rates. 
This is somewhat surprising at first, since it seems to us 
no previous work has explicitly noticed this point [39]. A 
plausible explanation is as follows: While e-h scattering 
is important for electrons and holes to effectively move 
together and thus maintain the validity of the ambipolar 
diffusion approximation, e-h scattering itself should not 
affect the diffusivity of the e-h ensemble since such scat- 
tering drives only internal dynamics. Rather, scatterings 
of the e-h ensemble with the ambient determines its dif- 
fusive capability in the ambient, which is the LO phonon 
subsystem in our model. 

To close this section we point out that the results de- 
rived for the limiting case of a predominant e-h scattering 
in Subsection V A and with standard ADA procedure in 
this subsection agree with each other. First of all, the 
derived diffusion equations for plasma density and tem- 
perature are identical. At the same time two approaches 
produce the same expressions for the DCs. The proof is 
easily seen by comparing Eq. (47) with its counterpart 
Eqs. (D4-D5) in Appendix D where explicit expressions 
for all the DCs are available for the single-component 
case. The results for the ambipolar temperature current 
can be derived in exactly the same way as in Section IV 
and thus will not be shown. Finally, we mention that 
we refer to the DCs for the single-component case simply 
as the ambipolar diffusion coefficients without distinction 
afterwards. 


VI. CARRIER GENERATION AND 
RECOMBINATION: 

CLOSURE OF fc-RESOLVED POLARIZATION 
HIERARCHY 


term and the spontaneous recombination term. First, 
we neglect the detailed carrier capture kinetics that car- 
riers undergo when entering the active region from the 
electrodes. Instead, an empirical model is adopted by 
assuming instantaneous carrier capture process for the 
injection current density J with quantum efficiency rj N 
and t]e . The factor 0 represents the spatial profile of the 
pumping current and A E g is the bandgap offset between 
the quantum well and the barrier material. Second, the 
spontaneous recombination term assumes a bilinear form 
as we ignore the correlation between electrons and holes. 
The stimulated interaction terms in the density and en- 
ergy equations are given by 

2 ^ 

R N\stim — R y Rstimjky r ) > (^I) 

k 

Rwlstim ~ R ^ 2 m Rstim(k,r) , (52) 

k 

where R s tim{k,r) is given by the last term in Eq. (5), 
which contains p(fc,r) explicitly. 

To close this hierarchy, approximation to the p(k,r ) 
is necessary. This issue has been addressed in the past 
for different special cases in two w^ays. The first is to 
eliminate the polarization equation adiabatically as done 
in Ref. [16] for the total density equation and in Ref. 
[40] for the kinetic energy equation. An alternative ap- 
proach is to replace the sum of the fc-resolved polarization 
with an effective polarization [11]. We will outline both 
approaches with slight generalization beyond what has 
been published in either case. 

We begin with the adiabatic elimination of the polar- 
ization. As was done in Ref. [41], the same adiabatic elim- 
ination can be performed in the presence of non-diagonal 
scattering terms in the p(fc,r) equation. To proceed, we 
introduce the slowly varying temporal amplitudes for the 
laser field and polarization: 


£{r, t) = i [E(r , t) + c.c] , (53) 

p{k,r) = | [ p(k,r ) e~ iUot + c.c.] . (54) 

Then the jfc-resolved polarization equation is rewTitten 
under the rotating-wave approximations as 


As mentioned, the moment equations still depend on 
k - resolved polarization p(fc,r) through the stimulated 
interaction contributions to the density equa- 

tion (41) and Rw\stim to the energy equation (43). The 
complete source contributions are specified as follows: 


d t p{k,r) = - [r 0 (fc) + iS k ]p(k,r) 

-i^-V{k)E 

a 




r t {k,k')- -v{h)v t>k -kf 


(55) 

P(k',r) , 


r n = 2 ^® - lnr N - B sp N 2 + R N \ Mm , (49) 

Rw = ^ AE g - t nr W - B sp NW + RwUtim ,(50) 

where a couple of approximations are made to obtain 
the analytical form of expressions for both the injection 


where 

M fc = E g {T L ) + h 2 k 2 /(2m r ) + Se e + Se k - ,(56) 

V{k) = r k + ti-l- (57) 
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After introducing the fc-resolved susceptibility Xk by def- 
inition: 


p(Kr)=soe b V^E, (58) 

Vk 

we can write the polarization equation under the adia- 
batic elimination approximation in the form: 


Xk = X° k + Y,W,k') Xk . , (59) 

fc' 

where the kernel JC and the zerot.h-order solution are 
given by 


xl 


K,(k,k') 


l^l 2 P(fc) 

* [r o (fc) + 

fir i{k t k’) - iV tM -wV(k) 
h[T 0 (k) + i&k] 


(60) 

(61) 


approach with the energy equation was partly addressed 
in Ref. [17] for the free-carrier case with a phenomenolog- 
ical polarization decay constant. In the following, we will 
follow the EBE approach to construct the kinetic energy 
equation. 

The EBE approach retains an effective (total) polar- 
ization equation by microscopically computing the total 
susceptibility and reconstructing the total polarization 
P(r,£). As a result, an equation for P(r,t) is obtained 
which is similar in form to the standard non-diagonal 
Bloch equation for a two-level system. The procedure 
has been described in detail in Ref. [11] and will not 
be repeated here. As can be easily seen, the total den- 
sity equation depends only on the total polarization after 
summing over h [11]. Therefore, P/v| stim will no longer 
depend on p(fc,r). The remaining difficulty is to deal 
with the total kinetic energy equations which still depend 
on fc-resolved polarization. 

Using the definition P a = (2/V) p* k p(k , r), we ob- 

tain from Eq. (55): 


This equation can be solved either by matrix inversion or 
by an approximation [41] similar to the Pade method as 
in the case of relaxation rate approximation [16]. Then 
the final expressions for the stimulated interaction con- 
tributions to density and energy changes can be written 
as 


fc 

(64) 

k 


^A r | stim 
Pw\stim 


k 

2fi 


(62) 


£n£t fi 2 (A: 2 ~ k' ) X*k v Xk 
v.fc-fc' 


k,k' 


2 m r pk 


Pk> 




2 h 2m r 

k 


\E\ 2 • (63) 


While adiabatic elimination of the polarization leads to 
a simple closure of the total set of equations, the resul- 
tant equations have a severe deficiency, especially in the 
presence of the spatial inhomogeneity. As was discussed 
in detail in Ref. [11], the reason is that the adiabatic 
elimination completely neglects gain dispersion. This the 
reason that some alternative time-dependent polarization 
equation was sought. One remedy is the so-called Effec- 
tive Bloch equation (EBE) approach as discussed in de- 
tail in Ref. [11]. The issue of how to combine the EBE 


where S 0 = 5 fc I*.=o - v s,k-k’- From Eq. (64), we can 
construct the following relationship: 


E'dtPa + Ed t P* a = iSo [EPa* - E* P a ] 

- * E ll^ E *P(k, r) - VkEp*(k, r )} .(65) 

TlV 2 Tfl f 

k 

Using this relation, we have 


Rw I stim — -Re{ E*(iSoP a + dtP Q } + Rw,ca ■ (66) 

While the first term now depends only on total polar- 
ization and its derivative, the last term, standing for 
the Coulomb-assisted (CA) energy change, unfortunately 
still depends on p(fc,r): 


J 


a h 2 k 2 

Rw,CA = r^E^ fc - fc 'T^ _ ^ fe, ’ r)p * (fc ’ r)_P * (fe,,r ^ (fc ’ r ^ 


2fiY 
i 

2 W 


k,k f 


Y, — P(k',T )V,, k - k ,p-(k,r) 


k,k’ 


2m r 


2m. 
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At this stage, we use the adiabatic elimination results for 
p(fc,r) to close the set of equations, since gain disper- 
sion has now been included with the effective polarization 
equations. 

Expressing the right hand side of the total energy equa- 
tion in terms of d t P was first done in Ref. [17]. Com- 
paring the first term of Eq. (66) with the correspond- 
ing equation in Ref. [17], we find that the term propor- 
tional to the polarization decay rate is absent in our ex- 
pression. This is because there the relaxation rate is a 
^-independent constant introduced phenomenologically. 
Here we have taken into account all (both diagonal and 
non-diagonal) scattering terms in the polarization equa- 
tions. The sum rule for the total scattering terms leads to 
the disappearance of the linear decay term that survived 
the ^-summation in Ref. [17] by contrast. 

Finally, plasma heating due to stimulated interactions 
results in a corresponding temperature change which is 
given by 


dfT\stim — JW T R\V | stim jN^Nlstim * (^^0 

To summarize this subsection, we have outlined two 
approaches to close the set of equations for the total ki- 
netic (thermal) energy and carrier density by approxi- 
mating the polarization equations in two different ways. 


VII. LASER FIELD EQUATION AND LATTICE 
TEMPERATURE EQUATION 

In this section, we specify Maxwell’s equation in a more 
concrete form. Our derivation will include frequency de- 
pendence of the background refractive index and the so- 
called thermal lensing effect, where the background index 
depends on lattice temperature. After making the slowly- 
varying envelop approximation to Eq. (1), it is written 
as 


— V’E— if E - = ^-JP. + A) + — fan , 

c 2 C 2 SoC 2 €o c z 

( 68 ) 

where the envelope functions P a and P& are defined 
through 


Wn 


2iuo , 


v a , b = 1 [P a , b e- iuot + c.c.] . (69) 

We assume that the background polarization satisfies 
the following “constitutive relation” in frequency space 
(with tilde added to the top of a variable): 


Pfe(w) = e 0 Xb ( lj,T l )E(u) , (70) 

where Tl stands for lattice temperature and Xb is in gen- 
eral complex with real and imaginary parts defined as 


usual: Xb = x' b + l X b ■ We Taylor-expand Xb around a 
given lattice temperature and the reference frequency 
uj 0 : 


Xb(^T L ) = x° b + (T-T° L )^ + (u J -uo)^+--. , (71) 

where = Xi»(wo,r°) = M 2 -l]+tX°" and we assume 
that there is no temperature and frequency dependences 
of the imaginary part of Xb (absorption). Substituting 
the Taylor-expansion into Eq. (70) and Fourier-transform 
the resultant expression to time domain, we will obtain 
the time-domain relation between P b (t) and E(t). After 
straightforward algebraic manipulation, the equation for 
the slowly varying envelope is written as follows: 


2K 


■(V‘ + K~)1 1 




A 


= ^Pa-ZrE+^-lkE 172) 

2 CqC 2 2 2 


where we have adopted the following shorthand notation: 


K 

8t 


Oib 


u) 0 nl 

3 

c 

1 . ® TLb 

2 (T l - Tl) dn b 
n° dT h ’ 


K 


o" 

xj 

K) 2 


(73a) 

(73b) 

(73c) 

(73d) 


The physical meanings of these parameters are obvi- 
ous: Pt describes the index change with temperature 
and accounts for the thermal lensing effects observed in 
high power laser operation. describes the background 
index dispersion and thus ,n° gives the group velocity 
index. Lastly, a b describes the background absorption. 

To complete our hydrodynamic model, we include the 
lattice temperature equation, which is modified from Ref. 

[4]: 


d t T L - d r {K ph C£d r T L ) = -y a {T L -T a ) 

+ C~ft (dtW\ L o + twojnrN + . (74) 

This equation takes into account equilibration to the am- 
bient temperature T a with a phenomenological rate 7 a , 
heat transfer from the EHP to the lattice due to c-LO 
scattering, temperature rise due to energy transfer from 
nonradiative recombination of e-h pairs, and certainly 
Joule heating by current injection into an active device 
wffth a congregate resistance R . Only phonon part of 
the specific heat C p h of the material is responsible for 
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lattice temperature change, and the negligibly small dif- 
ference between constant-pressure and constant-volume 
specific heat of the semiconductor material is ignored. In 
addition, a heat conduction term due to phonon diffu- 
sion is added with a thermal conductivity K p k . Finally, 
Eq. (74), together with Eqs. (41), (43), (72), the effec- 
tive Bloch equations (not listed here, but see [11]), and 
appropriate boundary conditions forms a complete hy- 
drodynamic laser theory. This set of equations describes 
couplings between lasing and heating self-consistently for 
a spatially inhomogeneous semiconductor laser. 


VIII. CONCLUDING REMARKS 

In conclusion, we have derived a coupled diffusion 
model (CDM) for the hydrodynamic variables — carrier 
densities and temperatures — for the electron-hole plasma 
in a semiconductor quantum well optical device. The 
derivation is self-consistent in the sense that the CDM 
is coupled to the optical polarization and the laser field, 
such that all relevant processes and variables are treated 
on the same footing. The major difference of this article 
from the earlier work is that carrier-carrier and carrier- 
LO phonon scatterings are treated explicitly in the mo- 
ment equations when deriving the energy and momentum 
decay rates. As such, all diffusion coefficients are given 
as functions of the densities and temperatures via mi- 
croscopical contributions. There are a few consequences 
wdth such a derivation. First, the resultant CDM not only 
consists of self-diffusion terms as most phenomenological 
models do, but also contains mutual-diffusion terms be- 
tween electron variables and hole variables and between 
temperatures and densities. Second, many-body effects 
of Hartree-Fock type appear in all of the diffusion coef- 
ficients very naturally, leading to modification of the co- 
efficients, which wall be discussed in connection with the 
numerical results for the coefficients in a subsequent ar- 
ticle. Furthermore, our explicit treatment leads to some 
new perspectives about the single-component reduction 
and the well-known ambipolar diffusion approximation. 
Finally, such a detailed treatment leads to the conclusion 
that the ambipolar DCs do not dependent on e-h scatter- 
ing. This work enables us to obtain and analyze all the 
DCs for an optical device in the general two-component 
case and under the ambipolar diffusion approximation 
from a microscopic point of viewc Such detailed numeri- 
cal results will be presented in a subsequent article. 

Another insight we gain from this study is the role of 
electron-hole correlation. The degree of such correlation 
determines whether the plasma subsystem can be ade- 
quately described as an ensemble of correlated pairs or a 
gaseous state of individual species, wffiich in turn requires 
either an effective single-component or two-component 
description. Within the scope of this work there are two 
distinct ways to establish such correlation: one is by e-h 
scattering in k space, and the other is by static Coulomb 
interaction in real space. These two types of correlation 


are responsible for the two types of single-component 
reduction, respectively. The first type of reduction de- 
scribed in Subsec. V A relies solely on &-space correlation 
(i.e.y e-h scattering), and therefore the internal field that 
leads to real-space e-h correlation is not needed. For the 
second type of reduction shown in Subsec. VB, real-space 
correlation is sufficient, and no assumption for e-h scat- 
tering and c-LO scattering is necessary. In reality, both 
real-space and A;-space correlations exist and are responsi- 
ble for the establishment of an effective single-component 
behavior in most cases. We note that the ambipolar diffu- 
sion approximation has not been derived self-consistently 
in the literature so far. Theoretical treatments thus far 
are mere adaptation to experimental observations, rather 
than rigorous derivation. The combined usage of these 
two types of correlation may eventually lead to a system- 
atic derivation of the ambipolar diffusion approximation. 
This will help elucidate when the ambipolar diffusion ap- 
proximation is out of question. 

Finally, a few' words about the plasma model (ignoring 
the excitons) are in order. As mentioned at the begin- 
ning of Section II, a plausible argument that supports the 
plasma model is the relatively high density in a typical 
laser. While this certainly holds in the phenomenological 
rate equation model which lacks of spatial dependence, 
this becomes a much weaker argument in a spatially in- 
homogeneous laser. It is w^ell knovm that carrier density 
is not uniform in the lateral direction of an edge-emitting 
laser and in the transverse direction of a surface-emitting 
laser due to pumping and carrier diffusion process and, 
to a lesser degree, due to the interaction wdth laser light 
field. While at the center of such a nonuniform distribu- 
tion, the density is certainly above the Mott density, but 
the density decreases to zero as we move towards the de- 
vice edges. Typically at about the half-value position of 
the distribution, the density drops below r the Mott den- 
sity. Such low density regions cause two problems: First, 
bandgap renormalization at high density level leads to 
the spectral overlap of the gain peak, 'where laser is reg- 
ularly designed to operate, wfith the exciton peak of the 
low r density region. Thus the laser field incurs strong ex- 
citonic absorption. This is especiall}' true for gain guided 
devices where the field tends to spread more into low den- 
sity regions. For a very short-cavity device wffiere mode 
spacing is quite large, w r e may be able to design a laser to 
avoid s\ich overlap of the laser cavity modes with the ex- 
citonic features. This becomes difficult for regular edge- 
or surface-emitting lasers where the longitudinal or trans- 
verse modes become closely spaced such that some modes 
wall always fall near the excitonic features. The effect of 
such absorption is expected to be weak in a strongly in- 
dex guided device, nevertheless. The second problem is 
related to the treatment of momentum and energy de- 
cay processes wffiich affect carrier diffusion. The plasma 
model is no longer valid in the low' density regions, as 
noted. Since the diffusion processes in such regions play 
a role for the overall carrier distribution and in particular 
for the density at the center, an appropriate treatment 
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of the e-h subsystem is thus necessary. Such a treatment 
would require a theory that takes into account all the 
intermediate situations from an excitonic solid, or con- 
densate, to a pure plasma phase. This is obviously an 
issue beyond the scope of this work, but one needs to be 
aware of this issue when dealing with a strongly inhomo- 
geneous lasing system. 


rates. Similar rates for the Maxwell distributions have 
been obtained in Ref. [37]. 

We assume that each band is in quasiequilibrium and 
described by a drifted Fermi-Dirac distribution function, 
n% = /fcLfco , as given in Eq. (17), where a = e,/i. Tsing 
the second Born approximation [44], the change rate in 
the distributions due to e-h scattering can be written as, 
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a«zu=f|V E (ad 

where V s , q is the Fourier component of the screened 
Coulomb potential, in the delta function, which 

stands for the difference in the total kinetic energy be- 
tween in- and out-scattering, is given by e% + e £ <+ q - £& ~ 
e k+q where a ^ /?, and T is defined as 


APPENDIX A: MOMENTUM AND ENERGY 
RELAXATION RATES DUE TO e-h 
SCATTERING 

It is well known that ultrafast e-e and h-h scatterings 
lead to carrier thermalization with a rate around 50 fem- 
tosecond at typical lasing density [20, 22, 24, 27, 43]. 
However, the role and consequences of e-h scattering is 
rather uncfear. In this appendix, we will study such a 
consequence of the e-h scattering within the second Born 
approximation. Specifically, we will show that the DB 
requirement of the microscopic e-h scattering processes 
lead to equal temperatures and drift velocities for elec- 
trons and holes. Linearization around the DB state al- 
lows us to obtain the momentum and energy relaxation 


T = K + y k ,(l-n%)(l-n 0 kl+g ) 

- n k n i+q( l ~ n fc')(! “ n k+q ) ■ ( A2 ) 

If the EHP is in quasiequilibrium, the DB condition 
requires that T = 0. Using 1 — = n£ exp[(e£_ fc a — 

H a F )/k B T% we can factorize T into the following form 




, 3 




where Q = n% +q n 0 ,(l-n%)(l-n%, +<} ) and u a = hk%lm a 
is the drift velocity. Using energy conservation, the above 
expression can be rearranged as 


T 


= e { 


1 — exp 


k + q 


eg ~ fa? /_1_ 

k B \T a 



(A3) 


by adding and subtracting the term hq ■ u a jk B T 3 . Since 
Equation (A2) is valid for arbitrary k and q, and Q con- 
tains no algebraic roots, it is then straightforward to con- 
clude that detailed balance leads to 

T e = T h , (A4) 

u e = u h . (A5) 

These equalities are the result of the DB requirement un- 
der the condition of strong e-h scattering and have been 
used for the single-component reduction in that limit- 
ing case. They are duplicated in Subsection V A as Eqs. 
(38-39). 

Now we consider the deviation of the EHP from the 


DB state. We make a linear Taylor expansion for the 
exponential term inside the curly' bracket in T , and the 
resultant expression is given below: 


r w -Gx ( T a - t 0 ) - g 2 ■ ( u a - u 3 ) . 


Plugging this linearized expression into Eq. (Al), we thus 
obtain the momentum and energy relaxation rates. 

The momentum relaxation rate due to e-h scattering, 
as appeared in Eq. (14), is defined as 
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d t P a \ eh = lThkd t nZ\eh (A6) 

k 

= (l) s E Y v --- S(A ‘“‘ )T ’ >k 

' ' fc,fc',q^O 

= -(§)’ E 

' ' fc,fc',g^O 

x [£, (r“ - T 0 ) + s 2 • (u“ - «*)] . 

We notice that there are two terms in the momentum 
relaxation, but ensuing proof shows that the coefficient 
of the temperature difference term ( T a — T 3 ) vanishes 
under the DB condition. This leads to the equation for 
momentum relaxation due to e-h scattering given in Eq. 
(29). To prove this, we notice that the prefactor of the 
term with ( T a — T 3 ) is proportional to 


£ V 2 q S (Ae eh ) n° k+q n 0 ,( 1 - n£)(l - n 0 , +q ) 

fc,fc' r g/0 

-e%-hq u a u _ 

k B 

£ V s 2 ,/ (Ae e ft) tf + X( 1 - /£)(! - /£+,) 

k,k’ ,q ^ 0 


X -~~r (k + k%) , 


k]3 


(A7) 


which are obtained by using Eq. (A5) inside the delta 
function after making the following translational trans- 
formations: 


fc — k jT) — k \ , 

fe' - k 0 D = k 2 , 

then change the notation back, ie., Aq — » A: and A: 2 A:'. 

Due to inversion symmetry in the prefactor of (k + fcg), 
only the term proportional to kp in Eq. (A7) survives. 
We rewrite its prefactor (xks) below: 

£ V 2 q S (Ae eh ) j£(l-/g)(l-/£ +fl ) (e£+, ~ eg) 

k t k J ,q^0 

We recognize that the expression denotes the inte- 
grated energy difference between in-scattering and out- 
scattering with respect to states \ak) and \j3k r + q) un- 
der detailed balance, and thus vanishes. Therefore, the 
momentum relaxation term is reduced to a simple rate 
equation form: 


7 eh 




k,k',qjt 0 


27 r 

h 


V s 2 q 6(ACek) 


(A9) 


hq x fokx 




This can be obtained by making the same replacements 
and using the inversion symmetry argument. We also as- 
sumed that the direction of u a - u 3 is taken as the x 
direction. The positive-definiteness of the above relax- 
ation rate can be verified easily after accounting for the 
isotropy of the EHP. Then we may replace q x k x in the 
summation by k q/ 2, which is related to (e£ +<? -cg)-cj. 
Since the energy difference term vanishes after integra- 
tion, the momentum relaxation rate takes the neat form: 


7e/i = 


2 m r k 


^-( 1 )’ y ■ 

tP _ rP 


^V 2 q 5(Ae eh ) 


x/ fc %/^(i-/ fc Q )(i-/ fc % 9 K- 


(A10) 


Next, we consider the energy relaxation term appeared 
in Eq. (15) which is specifically given as 


d t E° 


I eh — 


2 E^ a "*i«‘ 


S “ 2m a 

= (I)’ E 

^ 1*' f T- 


k,k’,q?0 
3 


27 r 


ft 5 ’ 9 


V 2 <5(Ae e/l )T 


--(I) E 


c e/i) 


x[^i (r tt -r*) + 0 2 •(«*-«*)] 


(All) 

n 2 k 2 

2m a 

h 2 k 2 
2m a 


Through exactly the same manipulation and arguments 
as we did for the momentum relaxation, the energy re- 
laxation term cab be written in the following form, 


9tE° | eh = -T? h {T a -T 0 ) 

-m r % h u a ■ (u a - u B ) , (A12) 

where the energy relaxation rate due to temperature dif- 
ference is given by 


r 


Of 

eh 


(I)’ E fo<^> 

^ ' k,k f ,q^0 

*/*%/*' (* - /fe )( J - 


(A13) 


a _ f oc 
e k 

k B T*TP ' 


d t P a \eh = -m r j e h (tx a -u 3 ) . (A8) Before closing this appendix, it is worth making two 

observations: (1) the relaxation rate equation for energy, 
Eq. (A12), is consistent with our previous results, Eq. 


where the momentum relaxation rate is given by 
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(27), and (A8). Thus, we prove, within linear expansion 
treatment, that 


dtW a left = -r“ ft {T a - T 0 ) . (A14) 

The corresponding temperature change rate due to e : 
h scattering can be written, according to Eq. (34), as 
follow's: 


d t T a U = j^d t W a \ eh = -7r« ( T a - T 0 ) , (A15) 

w'here the relaxation rate is given by 7j& = jw^eh- Fur " 
thermore, we point out that the relaxation rates are dif- 
ferent for electrons and holes. The reason is twofold: first 
and foremost, the T® h factor is inversely proportional to 
the carrier mass. As compared to Eq. (A10), we immedi- 
ately see a factor m T jm a is recovered for Eq. (A13), given 
that 7 e h is carrier mass independent. Second, the j (y fac- 
tor depends on the carrier mass in the quantum regime. 
Therefore, the temperature relaxation rate is inversely 
proportional to the carrier mass at low density. (2) the 
second equality, Eq. (A5), is the same as what is required 
in the ambipolar diffusion approximation [15, 16]. The 
requirement is consistent with the intuitive understand- 
ing of incoherent collisions between different species in 
the EHP. These interactions mean additional frictional 
force between the oppositely charged species, in addition 
to the force relative to the ambient, caused by interacting 
with mainly LO phonons. Such frictional drag may lead 
to a reversal in the drift direction of carriers, as corrob- 
orated in the study of the negative mobility of minority 
carriers in semiconductor quantum wells [29-34]. By con- 
trast, the ambipolar diffusion approximation is deduced 
by use of charge neutrality, which is a static condition. 
Interestingly, given a neutral initial condition for a sys- 
tem, it will evolve and preserve charge neutrality since 
there is no charge separation as the oppositely charged 
components move at the same velocity, as governed by 
the equation of continuity, that is Eq. (41), and Eq. (A5). 
Therefore, having ultrafast carrier scattering in the EHP, 
the ambipolar diffusion approximation is redundant for 
an initially neutral system. 


APPENDIX B: MOMENTUM AND ENERGY 
RELAXATION RATES DUE TO c-LO 
SCATTERING 

Interaction of the EHP with the host semiconductor 
crystal is dominated by the inelastic carrier-LO phonon 
scattering at room temperature, while collisions with 
other phonon branches are relatively weak and neglected. 
In semiconductor quantum well structures, as the w T ell 
width is reduced, phonon modes could become confined 
in the growth direction and interface modes are intro- 
duced for small enough width [45]. On the other hand, 
if the quantum well is not too narrow, calculation of the 


scattering rates using bulk modes produces similar re- 
sults to that obtained by incorporating both the confined 
and the interface modes for a semiconductor quantum 
well, if no mode-specific physics is concerned [46, 47]. In 
this work, we therefore use bulk LO phonon modes of 
the Einstein model for the 8 nm quantum well structure. 
Under the assumption of quasiequilibrium for the EHP 
and equilibrium for phonons, the momentum and energy 
relaxation rates due to LO phonon scattering, d t P a \ L o 
and d t W a \io are worked out in this appendix. 

First, let us put forth the necessary microscopic ingre- 
dients [45]. The rate of change in carrier distribution due 
to c-LO phonon scattering, according to Fermi’s golden 
rule, is given by 


dtnl\LO = 2 [<*( A <LK+,(l-n2) 

<h<lz 

-6{ Ae^o) ng(l-ng_„)] 

+ y£|tfa| 2 [«(A^ 0 K-,(l-"2) 

-<5(A e + 0 )ng(l-ng + „)] , (Bl) 

where phonon wavevector Q = (q, Q z ) is expressed in its 
in-plane component q and vertical component q z , lolo is 
the phonon circular frequency, H e and H a are the carrier- 
LO phonon scattering matrix elements for phonon emis- 
sion and absorption, Ac£ 0 = eg - eg +q + Nj lo , and 
A el Q = eg - eg - hui L o- Using the Frohlich Hamilto- 

nian, the matrix elements are given as 


\He\ 2 

\Ha\ 2 


e 2 hwLp 
2 eq£ p V 


I G{q z )\ 2 


Nq + 1 

Q 2 


e?j\ojLQ 

2eo e pE 


! G(q,)\ 


Q 7 


(B2a) 

(B2b) 


where l/e p = l/Soo ~ 1 fe s with e s (£«>) being the 
static (high-frequency) relative permittivity of the un- 
excited semiconductor. The overlapping integral, G{q z ), 
is defined as dzXi{~)Xj( z ) ex P(iQz z )> where Xi( z ) 
is the real and normalized carrier envelop wavefunc- 
tion in the i-th subband which has been always taken 
as the ground state in this work. The phonon den- 
sity is governed by Bose-Einstein distribution, Nq = 
1/ [exp {Plo^lo) - 1], with /3 L o = 1 AbTlo and Tio 
being the LO phonon temperature. Conservation of mo- 
mentum has already been explicitly considered in Eq. 
(Bl). 

Before calculating the momentum and energy relax- 
ation rates, we point out that we will still work in the 
linear regime, as in the preceding appendix, but now in 
the drift velocities themselves. (Given the femtosecond 
scattering time scale at typical lasing density [43, 48], it 
is reasonable to expect the drifted carrier distributions 
are small perturbations from the non-drift Fermi-Dirac 
distributions.) The linear expansion is well known as 
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(B3) 

where /'(eg) IS derivative of the Fermi-Dirac distri- 
bution function /(eg) with respect to the carrier energy. 

Now we are prepared to compute the rates. The mo- 
mentum relaxation term is defined in the same manner 
as in the carrier-carrier scattering case, that is 


Inversion symmetry leads to zero contribution from the 
non-drift part of the distribution, so we only need to 
focus on the linear term in k % in Eq. (B3). Lengthy but 
straightforward integration over the phonon wavevector 
Q produces a concise integrated solution. We only give 
a few intermediate steps in the deduction and the main 
result below: 


dtP^Lo^^ThkdtnZho- 

k 

I 


d t P a 


LO 


ne'uJLO 

(2n) 4 eoe p 


j ' dkhk J 


dq 


Fn(q) 


X (s(& ej 0 ) + /*)/£-,( 1 - /»_,) + + ')] /» U - /.)} 

- i(A 4o) {^£-(Nlo + - ft) + Kk t B T^ Wf - < A '“ + ] )1 A“V,U - /»"+.> 


where Nlo = Nq since the LO phonons are dispersion- 
less in the Einstein model and 


/ oo 

dq z \G{q z )\ 2 /{q 2 + ql) 

CO 

/ oo roo 

dzi / dz 2 Xi( z i)xi( z 2 ) 
-oo J — oo 


x exp(-g[zj - z 2 \)x\{ z 2)Xi{. z i) ■ ( B5 ) 


Then we evaluate the integrals in the momentum relax- 
ation term individually by using polar coordinates for 
both k and q. During the integration the direction of 
the drift wavevector is taken as in the x direction. Af- 
ter term-collecting and -canceling, the final momentum 
relaxation term is found to be 


d t P°\LO = -7 l 0 P a 


with the momentum relaxation rate given by 


(B6) 


Similar to the carrier scattering, it is found, after care- 
ful evaluation of the carrier energy relaxation due to LO 
phonon scattering, that the higher-order drift term can- 
cels out according to Eq. (28) as far as the carrier ther- 
mal energy is concerned. Thus we will only show the 
partial contribution from the non-drift first term in Eq. 
(B3) below. To assist computing the energy relaxation 
term, summation over the phonon wavevector in Eq. (Bl) 
can be carried out directly so that the procedure becomes 
more transparent. After tedious algebraic manipulations, 
the expression becomes 


dtn^LO = ci (I\ (k) {ATo/(cg - Hu LO ) [1 - f( e k)} 

-( N lo + 1 )/(c£) [1 - /(eg - hu LO )}} 

+ / 2 (A:){(A f LO + l)/(^ + ^Lo)[l-/(^)] 

-AW(eg) [1 - /(eg + toLo)]}) , (B9) 

where Ci = {m a e 2 u LO )/{ 2irti 2 £ 0 £ P ), while the two k- 
dependent integrals are given as 


pOO 

7lo = co J Io(k) {[Alo + /( e +)]/( e +) [1 - /( £ +)l 

+ [N L o + 1 - / W] /(«+) [1 - /(*+)]} dc , (B7) 

where Co = (m a e 2 u>Lo) / (8tt 2 h 2 £o£ p kBT) and e+ = e + 
X\u>to- We define k + = \/k 2 + 2m a ULo j h, which is used 
in the integral: 



F\ i(q) dq 

\J(2kqf - (q 2 + 2rn a u L o / hf 

F u (q) dq 

\J (2 ~kqf - ( q 2 - 2 m a Wio//i) 2 


/BIO) 

>(B11) 


rk + + k 

Io(k) = / 

Jk+-k 


q 2 F u (q) dq 


k \J (2/ctf) 2 - (? 2 - 2 maWLo/hj 2 


(B8) 


where k- = yj k 2 - 2m a wz,o/fi- The integrals are related 
by the transformation: I\ (k + ) = I 2 (k). Using the above 
equations for h t2 {k) and integrating over the polar angle 
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of wave vector fc, the carrier thermal energy relaxation is 
determined to be 

rOO 

d t W a \ LO = -ca / del 2 (k) 

x {(Nlo + 1 )/( e +) [1 _ /( c )] 

- N LO m [1 - /(£+)]} , (B12) 


component plasma in semiconductor quantum wells. Us- 
ing the solutions for electron ( a = e) and hole (a = h) 
momenta P a given in Eq. (30) , we can write the density 
current for the a-component in the form: 

J Na = - D N « x d T X - , (Cl) 

x q 


where c 2 — (m 2 a e 2 u 2 LO )/(2-K 2 h 3 £oe p ). Plugging the 

Fermi-Dirac distributions in, it can be shown that a fac- 
tor proportional to the temperature difference between 
the carriers and the host lattice exists. From equa- 
tion (34), -we can write the temperature change due to 
c-LO scattering in the form 

d t T a \ LO = j^d t W a \ LO 

= ~Tl 0 {T a - T lo ) , (B13) 

where ^lo is defined as follows: We note that the ex- 
pression inside the braces of integrand in (B12) can be 
rewritten as follows: 


where X € {N e , N h ,T e ,T h ) . We have introduced vari- 
ous DCs and conductivities which are listed as follows: 


Djv a N a ~ Ma 

jf 1 + Va)^N a + Tlfjc. 

D N « N& = Ma 

[(1 +T, a )Hfa+S*,' 

Djv<*T a — / i Q 

(1 + 7 1a )Sf.+H% a ] 

T) x a M a 

[(1 + r, a )H° 0 + S 0 e ] 

cr a = fi a e 

2 [N a (l + Va )-N 0 ] 


where a / /I and 


(C2) 

(C3) 

(C4) 

(C5) 

(C6) 


(Nlo + l)/(«+) [1 - /(f)] - AWW [1 - /(«+)] 

= (Nlo + l)/(c+) [1 - /(f)] 

x {1 - exp \((3 a ~ pLo)h^Lo)}} ■ (B14) 

Obviously the term is proportional to the temperature 
difference between the carriers and the phonons. There- 
fore, the constant can be explicitly written as 


r DO 

no = c 2 j^(N LO + 1) / dt I 2 (k)f(e+) [1 - /(f)] 

Jo 

(B15) 


exp [(fig - /?Lo)hu; LO )] 

T a - Tlo 


The constant is positively definitive for T a > Tlo - By 
taking the limit of T a -+ T LO in the above expression, 
we then obtain New+on’s cooling rate between electrons 
(or holes) and LO phonons: 


Tlo = c 2 j^(N L o + 

rOO 

x dt I 2 (k) f (t + ) [l — f (t)} . (B16) 

Here, the results are independent of the assumption that 
the lattice and LO phonons are in thermal equilibrium: 
Tl — T L o, despite that it is assumed so in this work. 


S% = d x W a + N a d x 6e Q , (X = N a ,T a ) , 

Hx = N a d x Se a , {X = . 

The factor /z a is defined in Eq. (31). The various rates 
used above have been given in the first two appendices. 
The corresponding temperature currents and associated 
DCs can be readily written down according to the rela- 
tion in Eq. (37) and will not be listed here. 


APPENDIX D: SELF- AND 
MUTUAL-DIFFUSION COEFFICIENTS UNDER 
THE SINGLE-COMPONENT APPROXIMATION 

In this appendix, all diffusion coefficients for the single 
component case are explicitly given in terms of material 
parameters, scattering rates, and distribution functions 
of the EHP in a semiconductor quantum well structure. 
We consider the limiting case of strong e-h scattering. 
This is true when the plasma density is high enough such 
that e-h scattering dominates over other scatterings for 
momentum relaxation. At the end, we introduce termi- 
nology for the diffusion coefficients, which shall be exten- 
sively used in the subsequent article. 

From the relations obtained in Subsection V A, it is 
easy to show that the currents for density and tempera- 
ture are given, respectively, by 


APPENDIX C: CARRIER DIFFUSION 
COEFFICIENTS IN GENERAL 
TWO-COMPONENT CASE 

In this appendix, we only give the expressions for 
density-related diffusion coefficients of a general two- 


jjv = Nu , (Dl) 

Jt = [2 jw (WIN) - jjv] Jn , (D2) 

where w'e have used Jt = j w J\v ~ js J.\ and Jw — 2u\\ 
with W being the total plasma thermal energy. 
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Recall that the drift velocity of plasma u is given by 
Eq. (40) and now we rewrite it in terms of the gradients 
of density and temperature as follows: 

u = ~^{d N W + Nd N 6e g )d r N 

-^(d T W + N&rSe g )d r T, (D3) 

where Se g = 5c e + 5e h is the total Hartree-Fock bandgap 
renormalization. As defined in Eqs. (44) and (45), the 
DCs are finally written as 


Dnn 

= n(d N W + Nd N Se g ) , 

(D4) 

Dnt 

= fi {d T W + Nd T 5e g ) , 

(D5) 

Dtn 

= {2iw{W/N) - j*a’] Dmn , 

(D6) 

Dtt 

= [2jw(W/N)-j N )D NT . 

(D7) 


In the two-dimensional case, analytical expressions for 
the energy and carrier density are given by 

W a = m ° ( ^ 2 ra) - r In (1 + e- 1 11 ) dx , (D8) 

J-0^% 

N ° = Tnak ?r In [1 + exp(/?a/i£)] , (D9) 

7T k Z 

The derivatives of W a with respect to N a and T a can 
also be written explicitly: 


d N «W a \ T <. = k B T a 6 a , (DIO) 

d T „W a \ N » = 2k B p a W a - k B N a 9 a , (Dll) 

where 0 a , defined as 

Oct — [1 + exp{— + ex P {PcxVf}] > (D12) 

can be considered as a degeneracy factor, which goes to 
one in the nondegenerate limit. Both (DIO) and (Dll) 
take their classical values in the nondegenerate limit. 

We also note from the definition of ^ [Eq. (31)] that 

/X = l/[m e 7£ 0 + mh'JLo] ( D13 ) 

when the c-LO scattering can be ignored compared with 
e-h scattering. In this case, all the diffusion coefficients 
defined in this appendix are the same as in Subsec- 
tion VB. This means that the two approaches lead to 
the same definition of ambipolar DCs. 

Before closing this appendix and for the benefit of dis- 
cussions in the text, we introduce the following terminol- 
ogy for the DCs. The coefficient is dubbed self-diffusion 
when it relates the gradient of a variable to its current, 
and mutual- diffusion otherwise. In addition, variable 
name is used to label the coefficient that relates to the 
variable gradient. For example, Dtn is called mutual- 
diffusion density coefficient accordingly. 
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